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On Diffusion Limited Deposition. 

A. Asselah, E.N.M. Cirillo, B. Scoppola, and E. Scoppola 


Abstract 

We propose a simple model of columnar growth through diffusion limited aggre¬ 
gation (DLA). Consider a graph Gat x N, where the basis has N vertices Gm '■= 
N}, and two vertices (x, h) and (x', h') are adjacent if \h — h'\ < 1. Consider 
there a simple random walk coming from infinity which deposits on a growing cluster 
as follows: the cluster is a collection of columns, and the height of the column first hit 
by the walk immediately grows by one unit. Thus, columns do not grow laterally. 

We prove that there is a critical time scale N/log{N) for the maximal height of 
the piles, i.e., there exist constants a < (3 such that the maximal pile height at time 
aN/\og{N) is of order log(A), while at time fiN/ \og{N) is larger than A^. This 
suggests that a monopolistic regime starts at such a time and only the highest pile 
goes on growing. If we rather consider a walk whose height-component goes down 
deterministically, the resulting ballistic deposition has maximal height of order log(A) 
at time N. 

These two deposition models, diffusive and ballistic, are also compared with uniform 
random allocation and Polya’s urn. 

AMS 2010 subject classifications; 60K35, 82B24, 60J45. 
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1 Introduction 


Motivation. A celebrated model of deposition via diffusion is proposed in the early 80 ’s 
by Witten and Sanders EH. The aggregate, denoted A(K), made of K sites of is built 
inductively as follows. Choose A(l) = {0} and assume A(K). Let dA(K) denote its outer 
boundary. Informally, launch a simple random walk, n h-)■ S(n), far away from the origin, 
and stop it when it reaches dA(K), say on random site V. We set A(K + 1) = A(K) U 
In other words, if tqa{k) is the time at which the walk hits dA{K), then for y G dA{K), 

P{A{K + 1) = A{K) U {y}\A{K)) = lim Px{S{toa{k)) = y\rdA{K) < oo). 


Simulations show that the cluster looks like a ramihed tree with long branches. Heuristically, 
the origin of reinforcement is clear. Think of the walk in terms of its radial component, which 
performs an almost symmetric one-dimensional walk, and its transverse component. Either 
the random walk sticks soon after reaching the outer radius of the cluster, and it has to 
settle on a tip, or it takes time before settling and its radial component diffuses, and has 
more chances to visit the extremal shells, hence increasing the probability of attaching a 
tip rather than an inside site. This explains reinforcement, but does not explain why this 
reinforcement is enough to produce a ramihed tree structure. It is clear also, at the heuristic 
level, that we face two problems: controlling the number of tips in the growing cluster, and 
controlling in a quantitative way the reinforcement of these tips. 

One natural way to measure the dimension of the cluster is to hnd the scaling of the 
radius of A{K), and look for d such that 

Radius(A(iL)) ~ (1.1) 


If A{K) were a ball, then d = d, and the conjecture is that d < d. Now, physicists have a 
much sharper conjecture 


In dimension 

Kesten in 
not too long. 


dr = d — 


d-l 

d+1 


( 1 , 2 ) 


2, dc = 5/3, and simulations give d = 1.7. 

da ESI HE] considers the problem, and shows that the arms of the cluster are 
More precisely, his result reads 


d> 


3/2 for d = 2 , 
2 for d = 3 , 
d/2 for d > 3, 


idr = 2-1/3) 
{dr = 5/2) 

{dr < d-3/5). 


(1.3) 


By reversing time, (see |2T] and assume d > 3) one writes the probability of adding Y = y 
to the cluster as 


P{A{K + 1) = A{K) U {y}\A{K)) 


f’y{'^dA{K) oo) 

'n,z^dA^^iffdA{K) = oo) 


(1.4) 


The difficulty is to have estimate on the escape probability when the set A is not a sphere, 
or some simple geometric shape. Let us mention an interesting result about holes in the 
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DLA cluster, where a hole is a hnite maximal connected subset of the complement of A(K). 
Erbez-Wagner [12] shows that in dimension two, almost surely the number of holes tends to 
inhnity with K. 

Barlow, Pemantle and Perkins in [7| study DLA on a regular d-ary tree where the con¬ 
ductance between edges joining generation n and n -|- 1 is for a < 1. These authors 
show that the inhnite cluster has a unique inhnite line of descent. Even though there is an 
explicit formula for the harmonic measure, the proof that r{A{K)) scales like K with normal 
fluctuations is non-trivial. 

Benjamini and Yadin in [5] propose another toy model for DLA. They consider a cylinder 
Gat X N, where the graph Gm has constant degree, N vertices, and is fast mixing: the mixing¬ 
time should be less than log^“'^(|GAr|) for some positive e (the class of d-regular random graphs 
works). They show that if we send H x |GAr| simple walks from inhnity, then the height of 
the aggregate is larger than id log(log(|Gjv|)) for any id and N large enough. 

There is a two-dimensional model, the Hastings-Levitov model, which takes advantage 
of the conformal invariance of two-dimensional brownian motion, and Riemman’s mapping 
Theorem to map the complement of the cluster into the complement of the unit disk, and 
then attach on the unit circle a stick at a random uniform angle. Recently Norris and Turner 
[26] have studied very precisely the limiting cluster obtained by iteration of randomly rotated 
conformal mappings. 

In a series of three recent papers, Amir, Angel, Benjamini and Kozma mw study DLA 
on Z with long-range random walks. The cluster is no longer connected, and they discover 
many phase transitions in the growth rate of the cluster according to the tail decay of the 
increment of the walk. 

Our model is a further simplihcation of Benjamini and Yadin’s model [9] in two ways: 
(i) no lateral hairs are produced, and (ii) the basis graph has no geometry. In our toy model 
of DLA, the radial component does a one-dimensional random walk, and the transverse 
component samples uniformly the section of our graph. Still we believe that our model is 
interesting, and one can answer some of the following questions in a quantitative way. 

• What is the origin of reinforcement? 

• What is the critical height to overcome ? 

• What are the different regimes in the cluster’s growth? 

Models. We shall consider two deposition models, diffusive deposition and ballistic depo¬ 
sition. 

We start with dehning diffusive deposition. Our graph is a half-cylinder Gat x N, where 
the basis has N vertices Gat := {1,..., N}, and two vertices (x, h) and {x', h') are adjacent 
if \h — h'\ = 1. The set Gat x {0} is called the ground. 

Let n I—>■ A{n) be the evolution of random subsets of Gat x N that we call the cluster. 
The cluster is built inductively with A(0) = x {0}. For an integer fc, the cluster A{k) is 
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made of columns, that is, 

N N 

A{k) = [J{ i} X {0,..., with = k. (1.5) 

i=\ i=l 

We shall write for simplicity A{k) = {ai{k ),..., aiq{k)). 

Assume that A{k) is built. We consider a simple random walk n Sn = (W„, 2^) on 
our graph. In other words, 

1. {Xn} an i.i.d. sequence uniformly distributed on Gn] 

2. {Zn+i — Zn} i.i.d. uniformly on { — 1,1}; 

3. the initial condition Zq is above the maximal height of the cluster A{k). For dehninet- 
ness we take Zq = max* (Jj(/c) + 1. 

The following rule of aggregation, or deposition, makes the cluster grow. The walk Sn, 
roams until it hits the cluster A{k). Let {X*, Z*) be the hitting site on A{k), and necessarily 
0 < Z* < ax*{k). We build A{k + 1) by increasing the height of column X* by one unit. 
That is 


ai{k + 1) = ai{k) for any i ^ X* and ax*{k + 1) = ax*{k) + 1 . 

We shall also say that the walk attaehes to the column, or pile, at X*. The walk with the 
aggregation rule is called an explorer. We shall denote by P the probability associated with 
this process. 

In diffusive deposition there are two relevant phenomena: one is diffusion, the other is 
deposition which happens instantly and this explains the name dijfusion limited deposition. 

Ballistic deposition is dehned similarly, with the same notation, but with a totally asym¬ 
metric walk {Zn+i — Zn = —!}• One could consider a continuum of biased models with a 
drift parameter. 

Definitions and notation. We use a,rj to denote conhgurations, i.e., o',r] & , a = 

(di, (72,..., (Ttv)- We also let |(t| ;= The symbol a will denote the conhguration 

obtained by ordering the components of a so that di > ^2 > ■ ■ ■ > dx- We call Ox the set 
of ordered conhgurations rj G , namely, such that rii>ri 2 >---> f]x- 

Given a conhguration a, we denote by ^(a) the height oceupation of a, i.e., 

N 

= • ( 1 - 6 ) 
i=l 

Note that = l<^l) and that C(o') = C(d). Given two conhgurations a and rj such 

that |(t| = |? 7 |, we say that a is more monopolistic than t], writing a t], when 

k k 

\/k = i,...,N 

i=l i=l 


(1,7) 
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Equivalently, one realizes rj from a by moving particles from the highest columns to the 
lowest ones. 

Urn models are paradigms of reinforcement phenomena (see for instance the survey m, 
and our deposition models actually can be stochastically compared with urns with N colors. 
We briefly recall Polya’s urn with N colors: starting with one ball of each color, at each 
unit time one draws a ball and put it back in the urn with an additional ball of the same 
color. Calling rji, with i = 1,..., N, the number of added balls of color i after \ri\ draws, the 
probability of drawing a ball of color i is 


Qf iv) = 


hi + 1 




( 1 . 8 ) 


We consider also a generalized urn by replacing the r.h.s. in (1.8) by /(hi)/ /(hj) where 


/ : N —)■ is a function such that /(O) = 1. When f{x) = + 1, we call the model the 

quadratic urn. When / = 1, we call the model the uniform random allocation and denote by 
{rj) = 1/N the corresponding probability of drawing a ball of color i at any time. Finally, 
we say that a process t h->■ a(t) is more monopolistic than process t i—)■ pit) if there is a 
coupling of the processes such that for any f > 0 we have a{t) >- pit), if this is the case 
initially. 


Main Results. In this Section we collect our main results. The hrst Theorem gives an 
estimate of the number of explorers necessary, in diffusive deposition, to form a cluster with 
at least one column proportional to a power of N. 


Theorem 1.1 Consider diffusive deposition. There are constants a < (3, such that almost 
surely, when N is large enough 


, aN 

maxcrd-—T— 
ieGjv ^log(iV) 


) <31og(iV), 


and there exists a positive constant x such that 


max (Tj 
i^G N 


f3N 

log(iV) 


) > N^- 


(1.9) 


( 1 . 10 ) 


In ballistic deposition, we prove that the growth of the height of the cluster is much slower. 
Indeed, it is unlikely that N explorers produce a column of height log(iV). 


Theorem 1.2 Consider ballistic deposition. There exists a positive constant A such that 
almost surely, when N is large enough 

max(Ji(A^) < 741og(A^). (l-H) 

i&GN 


Remark 1.3 For the radial component, we could have chosen — Z^] i.i.d. with some 

hnite range law without affecting our results. 
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Theorem ]_T occurs in a regime where less than N explorers are thrown in the graph. We 
call this the early regime which is to be thought of as the configurations where an additional 
explorer has good chances to settle on the ground x {0}. To motivate other results, let 
us explain the different steps leading to a column of height N^. A first step is to reach a 
subcritical height log(A^)/ log(log(A^)). We obtain that a large number of columns reach this 
height by comparison with random allocation. We obtain interesting comparison with other 
urns, with the observation that ballistic deposition looks like Polya’s urn with N colors, 
whereas diffusive deposition looks like a quadratic urn (see below ( |1.8 ) for the dehnition). 
Then, one of these subcritical columns reaches the critical height log(A^). Since our estimate 
requires the configuration to stay in the early regime, one has to bound the number of critical 
columns. We show that the number of critical columns is less than for some positive 

X, and this implies that the evolution remains in the early regime as long as the highest 
column has not crossed N^. We now can state our comparison result. 


Theorem 1.4 Both deposition models (diffusive and ballistic) are more monopolistic than 
Polya’s urn, which itself is more monopolistic than random allocation. 


The following corollary is a side result interesting on its own right which seems new, to the 
best of our knowledge. 

Proposition 1.5 Polya’s urn with N colors is monotone with respect to the order y. 


Related Models. There are many models of cluster growth similar in dehnition to DLA. 
They differ according to the law of Y, the site we add on the boundary of the cluster A. 
This can also be expressed according to the site, say X, from where the random walks are 
launched and lead to different phenomenology. 

• If X = 0, we rather dehne a dual model of erosion. The cluster represents the eroded 
materia, and A(0) = 0. Each new walk starts at 0, and settles on the hrst visited site 
outside the cluster (a site which we interpreted as being eroded). This is internal DLA, 
and was introduced by Meakin and Deutch in [23]. The cluster is spherical as was hrst 
seen Lawler, Bramson and Griheath in [22|. The huctuations were studied in in El E] 
and independently in [THl EHl 120] . 

• If X is uniformly drawn in the cluster, then Benjamini, Duminil-Copin, Kozma, Lucas 
in [5] show that the cluster is spherical. 

• If X is uniform on the boundary of the cluster, then this is the celebrated Eden model 
nn, which was proposed in the ’60, and studied hrst by Richardson [28] . 

• If particles do not erode immediately the materia, but do it with an exponential clock, 
and if they can be activated again when another walk stands on their site, this is 
Activated Random Walks. This model has been introduced by Spitzer in the 70, and 
much discussed in the physics literature as an example of self-organized criticality. 
This has been studied mathematically by Rolla and Sidoravicius [22] (and references 
therein), and recently by Sidoravicius and Teixera [30] among others. Recent ehorts 
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have focused on the case of an initial condition drawn from a product Poisson measure. 
As one tunes the density there is phase transition between settlement of explorers (in 
any hnite box), and their perpetual activity. 


Pictures and simulations. In order to illustrate our main results, we show some numer¬ 
ics. In particular, we emphasize the freezing phenomenon which leads to the monopolistic 
regime: after a given time the highets pile grows linearly catching all particles. We stress 
that simulations do not capture quantitative aspects of the problem (scalings or exponents), 
but serve merely as qualitative illustrations. 

For both the diffusive and the ballistic model we have simulated the systems for N = 
50,100, 200,..., 1000. For the diffusive model we have considered also the cases N = 
2000,4000,5000,6000,8000,10000. In all the cases averages have been computed over 10"^ 
independent realizations of the process. We have checked in all the cases that the sample is 
large enough! to get stable averages. 



4.. 8000 

7000 
.2 6000 
I 5000 
.3 4000 

8 3000 

^ 2000 
g, 1000 
^ 0 

0 2000 4000 6000 8000 10000 

graph size 



Figure 1: Numerical simulations for the diffusive model. Left panel: the highest column 
height is plotted as function of time (number of explorers). The hve plotted curves, from the 
left to the right, refer to = 100, 300, 500, 700, 900, respectively. Right panel: solid disks 
refer to the simulated highest column height at time N, namely, after N explorers have been 
sent, for different values of the size of the graph N. The solid line is an eye-guide obtained 
by plotting the htting function 0.498 x 


Simulations show clearly that the ballistic model reaches the monopolisitc regime much 
later than the diffusive one. Indeed in both models, compare the left panels in Figures and 
the height of the highest pile attains a linear behavior after an initial transient. This late 
time regime is the one in which all the particles are caught by the highest pile. Data show 
that the time length of the transient is much smaller in the diffusive model. 

Indeed, we 


We have also tested numerically our main results in Theorems 1.1 and 1.2 


have computed, by averaging over different realizations of the process, the typical height 
of the highest pile at time N. Rigorous results suggest that this quantity should scale as a 
power law in the diffusive case and logarithmically in the ballistic one. The related numerical 
results are shown in the right panels in Figures [T] and The qualitative agreement between 
simulations and theoretical results is striking. We stress again that the numerical results 
cannot be interpreted as a quantitave description of the model behavior since, for instance, 
too small values of the graph size N have been considered. 
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Figure 2: Numerical simulations for the ballistic model for N = 50,100, 200,..., 1000 with 
averages computed over 10^ independent realizations of the process. Left panel: the highest 
column height is plotted as function of time (number of explorers). The plotted curves, from 
the left to the right, refer to iV = 50,100, 200,..., 1000, respectively. Right panel: solid disks 
refer to the simulated highest column height at time N, namely, after N explorers have been 
sent, for different values of the size of the graph N. The solid line is an eye-guide obtained 
by plotting the function 1.957log(iV). 


Finally, in the diffusive case we have also tested numerically our results about the critical 
character of the time scale iV/log(A^). In Figure]^ we compare the highest column height 
measured at times N/ log(A^) and 2N/ log(iV). In the hrst case the numerical (solid circles) 
data can be perfectly htted by a logarithmic function. In the latter case, on the other hand, 
the poor logarithmic htting is opposed to a perfect power law one of the numerical data 
(solid squares). This result is in perfect agreement with the one proved in Theorem 0 and, 
in particular, it suggests that 1 < a < /3 < 2. We stress that our numerics cannot in any 
case be considered quantitative, indeed, we have no clue to state that, by considering larger 
sizes of the graph, our numerical results would be conhrmed. 



Figure 3: Numerical simulations for the diffusive model. Solid disks and squares refer, 
respectively, to the simulated highest column height at times iV/log(A^) and 2iV/log(A^) 
for different values of the size of the graph N. The solid line is an eye-guide obtained by 
plotting the htting function 1.847 x The two dotted lines are the graph of the two 

functions 0.724log(A^), 4.085log(iV), and 5.292log(A^). 
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Plan. The rest of the paper is organized as follows. In Section we present onr main tool, 
which is the probability an explorer hits the gronnd, as well as a henristic explanation for 
the logarithmic scale of the critical height. In Section we present comparison with urns 

is given in Section [T2 The proof of the Theorem 1.4 


models. The proof of Proposition 1.5 


is given in Section 3^ In Section]^ we establish our main tool, and related estimates. We 
study the very early regime in Section Then, we study the growth of cluster in Section 
and the reason why a large number of columns cannot overcome height log(A^). Finally, we 
gather all the needed estimates to prove Theorem IT in Section 


2 Key Tools and Sketch. 


The key to Theorem |1.1| is an estimate of the probability of attaching to a given column. 
Before, we need a lower bound on the probability of hitting hrst the ground. 


Lemma 2.1 Consider diffusive deposition with a configuration a such that |(j| < N/2 ex¬ 
plorers, 

1 iv II 

Pg{cr) := P( Explorer hits the ground | a) > exp + (l + 0(^) (2.1) 

1=1 


The time spent on the slab Gn x {0,..., ai} before touching the ground is typically af for a 
SRW, but only if the walk has good chances to cross the whole slab. Our key attachement 
estimate follows. 


Lemma 2.2 Consider diffusive deposition. Let a be a configuration such that |cr| < N/2. 
Then there exists a positive constant kd such that 


2 1 2 TV 

P( Explorer attaches pile i \ a) > —xexp (E‘"i(v+l))(l+0(^)). (2.2) 


1 = 1 


In a sense (2.2) and (2.1) are saying opposite things: the former inequality tells how easy it 
is to get trapped, whereas the latter tells how easy it is to reach the ground. 

Imagine a regime where ^ view of (2.1) the probability of hitting the 

ground would be small, and very likely the walk would not go below P, where P is such 
that 

5^ (a. - H)l ~ N. (2,3) 


In other words, P of (2.3) would play the role of an effective ground. We then replace (2.2) 
by the following estimate. 

Corollary 2.3 In the diffusive case, and for any positive H, 


N 


2 ^ 

P( Explorer attaches pile i\cr) > kd ——2:xexp + (2.4) 


X=1 
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Sketch. We wish to sketch heuristically the reason why log(iV) is the critical height at 
which a monopole forms. We hx a given column, say column 1, and we estimate the number 
of explorers needed to produce a given height. Lemma 2.2 allows us to bound this number by 
a sum of independent geometric variables, for which we know everything. Indeed, introduce 
Ti the number of explorers needed so that the height of our distinguished site reaches height 
1, that is Ti := inf{n > 0 : (Ji(n) = 1}. By induction, for any integer h knowing Th we 
dehne Th+i := inf{n > 0 : ai(r/i + n) — ai(r/i) = 1}. Assume now that we are in a regime 
where hitting the ground is likely. The estimate (2.2) means that for any integers h,n 

h\ 


P(Th+l > n\ Tl, . .. ,Ti) < (1 - 


(2.6) 


Let US now introduce independent geometric variables {r^, h > 1} with E[Th^]\ = 
Then, we will show that for any height H 


H 


H 




( 2 , 6 ) 


2 = 1 


2 = 1 


Recall that T < means that X explorers produce a column of height H at site 1. 

Now, we want to hnd X such that a given height H is likely to be reached. This would be 
the case if the probability that any distinguished site reaches height H is above 1/N. Thus, 
we look for X such that 


H 




2=1 


(2.7) 


Let us write X as N/f{N), and try to guess the size of f{N) which produces a monopole. 
Note also that for any H > f{N), 


H 


H 


N 


- E £) 6 < 


N 


i=f{N) 


i=f{N) 


KDf{N)' 


Thus, < A^//(A^)} imposes a constraint only on the hrst f{N) variables in (2.7). 

We then have to estimate f{N) such that 


fX) 


nY. 


2 = 1 


N - f{N) 


^ N 


We use now that {tj} are independent geometric variables 


N 


fW 


N 


fX) 


p{f, + --- + > n X. < 757m) n (l - (1 - ^) 


... ■ PW 

fW 

> (^1 - exp(-K£) 

2 = 1 

/ 1 \ fW 


2=1 


KdI \N/fiN) 


N 


PiN) 

ifmr- 


fim 

n 

2=1 


Kd 


P(N) 


( 2 . 8 ) 
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Now, using Stirling’s formula, we obtain 

/ 1 \ fW / \ 

P(ri H - h fff < X) > ^kd-^J = exp - log (—) X /(iV)j. (2.9) 

Thus, if f(N) is of order log(iV), it is likely that one monopole forms. 


3 Comparison with Urns 

In this section, we establish a coupling between our deposition processes and simpler ones 
which preserves a natural order on ordered conhgurations, to be dehned below. 

We consider growth evolntion on snch that at each unit time we add a unit height 
to a conhguration, say p, at a given site, say i, with a probability Piirj) which depends only 
on the value rji, and on the unordered set {rjjtj ^ i}. In this case, it is usefnl to reorder the 
indices through a permutation of the indices to obtain conhgurations whose heights are in 
decreasing order. We call p = {(pi(?7), • • • ,PAr(?7)), p G the law of the growth process. 


3.1 Comparing evolutions 

It will be important to compare conhgnrations with the same number of explorers. Our main 
results are the following. 

Proposition 3.1 Consider two processes t h-)■ rj{t) and t h-)■ a{t) on evolving, respec¬ 
tively, according to the laws p = (pi(-)) • • • )Pv(-)) and q = {qi{-), ■ ■ ■ iqNi,-)) ■ Assume that 
for any rj,a & Ojsi such that p ^ a we have that 

k k 

VA; = l,...,iV (3.1) 

i=\ i=\ 

Then, the process a{t) is more monopolistic than p{t),i.e. there is a coupling between the 
two processes such that a{t) >- pit) for any t. 


Lemma 3.2 Let {pi,. 
up to 1. Assume that 


Then, for any k = 1,.. 


.. ,Pn} and {qi,..., q^} two sets of positive numbers both summing 


Pi “ P2 “ ~ Pn' 

. ,N, we have that 


(3.2) 


k k 

> 5^Pn 

i=l i=l 


(3.3) 


Proof The proof is by induction on N. Assnme the Lemma is true with iV — 1 sets of 
positive numbers, and dehne the renormalized iV — 1 numbers 


1 


1 


Pi = Pi X 


1 -PAf’ 


and di = qi X 


f — qN 


(3.4) 
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The induction hypothesis states that for A; = 1 to — 1 

k k 

i=l i=l 


In other words, 


^ (1 — Qn) 

2=1 ^ ' 2=1 


Pi 


(3.5) 


The question is whether 1 — > 1 — Pat or Qn/pn < 1 which follows from (3.2) since 


1 


N 


N 


Qi ^ Qn 

i— > ypi — 

*=i tr 


Z—/ 1 


Qn 

Pn 


As a corollary of Proposition |3.1| and Lemma [3.2| we have the following result. 

Corollary 3.3 With the notation of Proposition [)0 , assume that for any r],a ^ On such 
that 7] -< a we have that 

Piih) / Qii^) 


< 


for i = 1,..., N — 1. 


Pi+iiv) 

Then, there is an order preserving coupling between ri{t) and cr{t). 


In order to prove Proposition |3.1| we need some notation and some simple observations. 
We dehne the action Aj : —)■ of adding one explorer to site j: {Ajri)i = rji + Sij. Note 

that Aj does not leave On invariant. 

Assume that rj E On and dehne 

I{ri) = {i E {2,...,N} : U {1}. (3.6) 

Also, for i G {1,... ,N}, let d{ri,i) = maxJ(p) n{l, ... ,i}. In other words, d{ri,i) is the last 
position of a height decrease up to position i. Note that for p G (Pa? we have 

Av = A{'n,i)V E On- (3.7) 

For rj E On, note that if i < j, then d{r],i) < d{r],j). Also, if i < j, then Ajp -< Aip. 

The main observation about ordering is the following. 


Lemma 3.4 Assume that p,o' E On with rj -< a. If i < j, then Ajp -< Aia. 


This lemma is based on the following simple observation. 


Lemma 3.5 Assume that r],^ E On with rj -< a. If some integers i < j satisfying d{ri,j) < 
i, are such that 

i i 

'^Pk = ^(Tk, (3.8) 

k=l k=l 


then i > d{a,j) 
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Proof of Lemma Assume for a moment that i > 1. Since rj -< a, we have 

2—1 2—1 2+1 2+1 


(Tfc, and Vk — ^ ^ ^k‘ 


(3.9) 


Rewrite now (3.8) as 




k=l 


i—l 


k=l 


k=l 


k=l 


i+1 i+1 

'^Vkff = CTi, and '^Vk- r]i+i = '^ak- Ui+i. 

k=l k=l k=l k=l 


Using (3.9), we have both that cxj < pi, and rji^i < (Tj+i. Since a G On, we have cxj+i < (Tj. 
Now, if i = 1, we have cxj = rji so that (jj < rji is again true. So that we reach 

Vi+i < < CTi < Vi- (3-10) 


Now d{r],j) < i < j means that rji = = ■■■ = Vji with (3.10) this implies that 

CTj = (Ti+i, and by induction, we reach that 


Tji (Ji C^2+l ■ ■ ■ O'j. 

These last equalities mean that i > d{a,j). 


(3.11) 


Proof of Lemma S.f. To simplify notation assume cr G On- We have already observed 
that for i < j, AjP -< Aip. Thus, we only need to prove that AjP -< AjCr. If d{ri,j) = j, 
then d{ri,j) > d{a,j), and the result is obvious. Assume henceforth that d{r],j) < j. If for 
all k = d{ri,j), ... ,j — 1, we have that 

k k 


2=1 


2=1 


then the result is also obvious. In the opposite case, let k in [d(?7, j), j[, be the hrst index for 
which we have 




2=1 


2=1 


then. Lemma [3.5| implies that k > d{a,j), and the lemma follows. I 

Proof of Proposition \ 3 . ij By way of induction, assume that up to time t, we have r]{t) -< cr{t). 
Draw a uniform random variable U in [0,1[, and dehne two random variables J, J* as follows: 

- if f/ G [piivii)) H - Pi-i{vit)),Piiv{i)) H- Piiv{i))[ then J = i {we set po = 0); 

- if U G [qi{(T{t)) H-h qj-i{(T{t)),qi{(T{t)) H-h qj{(T{t))[ then J* = j. 

Then, ( 3.1[ ) implies that J > J*. We set 

p{t + 1) = Ajp{t) and a{t + 1) = ^j*cr(t). (3.12) 


Then, Lemma 3.4 yields that p{t + 1) -< a{t + 1). 
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3.2 Comparing Polya’s Urn with Random Allocation 


By random allocation, we mean repeated draws of one out of N colors, labelled from 1 to 
N, uniformly at random. In other words, at each draw, the probability to pick up color i 
is 1/N. The law for Polya’s urn and random allocation are denoted respectively and 
with 


Va e 



+ 1 

N + Yli<N 


and 



Lemma 3.6 Polya’s urn with N colors is more monopolistic than random allocation of N 
colors. 


Proof. Note that for any 


which clearly holds since 


Proof of Proposition 1.5. 


<J,ri ^ On such that rj -< a 

(if (^) > (iffl) + l . 

qf+iW) - qY^.^{ri) ^*+1 + 1“ 

a G On- Thus, Corollary |3.3| implies the lemma. 
Note that ii p a and \p\ = |cr|, we have 


(3.13) 


VA; = 1,..., iV, 




+ Eti hi > + Eti 

^ + Ef=ih. + 


'^(ifi.(^)- 


This establishes Proposition 1.5 saying that Polya’s evolution with N colors preserves the 
order. I 


3.3 Comparing deposition models with Polya’s urn 


Recall the definition of ballistic and diffusive deposition given in Section [T} Denote their 
law, respectively, by p^ and p^. We show that both ballistic and diffusive deposition are 


more monopolistic than Polya’s urn, which is one of the statements of Theorem 1.4 
Assume for a moment the following lemma. 

Lemma 3.7 For any k = 1,... ,N and p G On 

k k k k 


(3.14) 


2 = 1 


2=1 


2=1 


2=1 


By Proposition |1.5[ Lemma 3T, and Proposition |3.1| we have that both ballistic and diffusive 
deposition are more monopolistic than Polya’s urn. 

To state a preliminary simple observation, we need more notation. For p G On, let 


pfk^iv) be the probability that the explorer hits site i at height k, for k eN. 
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Lemma 3.8 For any i, G {1,..., A^} and T]{i) > k > k' > 0 we have 

Pf,k'iv) < Pf,kih)^ ^fvU) > k pf^kiv) = Pf,kiv) 

Similarly, 

pSei’i) < */ ’Ki) > k p'J,(t)) = p“t())) 


(3.15) 


(3.16) 


Proof of Lemma 3_D_. We consider first the ballistic case. In view of Lemma 3^, we need to 
show that 

<1^ (h) 


Pf{h) > Vi + ^ 


pf+iiv) hm + 1 Qi+iiv) 


In order to prove (3.17), we need to show that 


> ESpf+i,.(h) 


hi +1 


hi+l + 1 


(3.17) 


(3.18) 


By Lemma 3.8, this ineqnality has the strnctnre 

ai + ■ ■ ■ + On ^ Ctm+l - 


> 


n 


n — m 


for n > m > 1 and ai > 02 > ■ ■ ■ > Un- The validity of snch an ineqnality is immediate once 
we let p = (om+i + ■ ■ ■ + On)/{n — m), note ai > ■ ■ ■ > Um > p, and write 


ai 


n 


ai + ■ ■ ■ + Om + p{n — m) pm + pin — m) 

- > -= p 


n 


n 


Finally, by nsing (3.17) and Lemma 3.2 the hrst of eqnations (3.14) follows immediately. 


The diffnsive case can be treated in the same way. This completes the proof of the lemma. 


Proof of Lemma 3^^ First we prove the lemma for the ballistic case. The lemma follows, 
since ^ 

pf,kiv) = explorer snrvives till k + 1) 


and P(the explorer snrvives till k) is an increasing fnnction of k. 

Now, we consider diffnsive deposition. First, assnme k — k' an even nnmber. To each 
path s' = {(a:)-, ^))}j=i,...,n hitting 77 in z at height k' in a time n we associate nniqnely a path 
s = {(a^j, ^j)}j=i,...,n, of the same length and, therefore, the same probability, hitting 77 in z 
at height k. Thns, the lemma follows since there exist other paths ending in z at height k. 
Given the path s', we constrnct s in the following way. Call zzi the time of last passage of s' 
throngh the intermediate height H = {k + k')/2. s is eqnal to s' np to time tzi while after 
7 Zi it nses opposite height increments than the original s', i.e. Zj+i — Zj = —— z'j) for 
all ni < j < n keeping the same horizontal increments. We therefore obtain a path ending 
in z and height k. Note that snch a path avoids 77 nntil it hits site z at height k, becanse 
77 is a nnion of colnmns. li k — k' is an odd nnmber we do a similar constrnction bnt we 
have to associate a set of paths s' of lenght zz to a single path s of lenght n — 1. The set 
is obtained considering together all the paths s' coinciding everywhere bnt the component 
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Xni+i, where ni is now the last hitting time of the level of the vertical process. Then 

the path s hitting rj on site i at height fc in a time n — 1 coicide with the paths s' up to 
time ni, say Xj = x'p Zj = z' for any j < ni and is specular to them after ni + 1 , that is 
Xj = Zj — Zj^i = —{z'j^^ — z'j) for any j = ni + 1,n — 1 so that Xn-i = i and this is 
the hrst hitting to p. Clearly the probability of s is larger than or equal to the sum of the 
probabilities of the paths s', since s is one step shorter, and the sum on is done only 

on - Cni+i(h) sites. 


4 Estimating Unit Growth 


In this section we discuss how heights grow. We consider the random walk Sn = (W„, ^n), 
and for an integer k, we call Hk the hrst time the walk reaches height k. In other words. 


Hk = inf{n > 0 : Zn = k}. (4.1) 

Since the X-component is uniform on the base, giving the conhguration a, the ordered one 
a, or the height occupation ( (dehned in ( 1 . 6 )) is equivalent, and we use Pg{<j) or Pg{C) 
indifferently to denote the probability an explorer hits the ground. Lemma |2.1| is obtained 
as a simple application of Jensen’s inequality whereas Lemma 2.2 requires Kesten-Kozlov- 
Spitzer representation of the local times d- 


Proof of Lemma |2.1[ For an integer k < ai let l{k) be the number of visits of height k 
by the random walk before Ho. We have the represenation 

PAO = (1 - |)'“’] = B[exp (f;!(A.)log (1 - |)). (4.2) 

k=l k=l 

Our hypothesis |(t| < N/2 implies that for k > 1, ^ Ci ^ -^/ 2 , and since log(l — x) > 

—X — for 0 < X < 1/2, we have using Jensen’s inequality 

P,(C)>exp(-f;£|;(4)](| + ^)). (4.3) 

k=l 

Now note that £'[/(/c)] = 2k. Indeed, the height of the random walk being a simple random 
walk on N, we have for A; < di, by conditionning on the hrst step 

E[l{k)\Zo = di] =E[l{k)\Zo = k] = l + ]^{E[l{k)\Zo = k + l] + E[l{k)\Zo = k - 1]) 

=1 + ^{E[l{k)\Zo = k]+ P{H, < Ho\Zo = k- l)E[l{k)\Zo = k]) (4.4) 

=1 + ]^{E[l{k)\Zo = k] + {l- \)E[l{k)\Zo = k]). 
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The equality E[l{k)] = 2k for k < ai follows at once. Note now that 


N 

'^2kCk = + 1). 

k>l i=l 


Finally, (2.1) follows as we note that 


Y,2k 

k>l 


- N 


k>l 



O 1*^1 + 1) 

N ^ N 
2=1 


On Kesten-Kozlov-Spitzer representation. Let u{h) be the number of up-crossings of 
height h before touching the base. In other words, we define m( 0) = 0 and for h > 0 


Ho 

u{h) = '^l{Zi.uZi)={h,h+i)- (4.5) 

2=1 

Similary, down-crossings of height h correspond to jumps from h to h — 1 before time Hq. 
One way to realize the random walk n i—>■ is to assign the sequence of up and down- 

crossings on each height. Thus, we consider i G N}, k E N} a. collection of i.i.d. 

geometric variables, with law P{^ = n) = 1/2'^^^ for n G N. Now, the sequence of up and 
down-crossings at height k is as follows: .^q up-crossings, then one down-crossing, then 
up-crossings, the one down-crossing, then ^2 up-crossings... and so on and so forth. The key 
observation is that each for i > 1, is preceded by an up-crossing of the height k — 1. In 
other words, 

li(/c— 1 ) 

=^o+ = ^0- 

i=i 

We set G{h) = k < h, z G N) the cx-held representing the choices of moves on the 

hrst h heights. Kesten-Kozlov-Spizter representation expresses the local times of Z in terms 
of the u. Thus, if l{k) represents the number of visits of height k before Hq, for a walk with 
starting level above di, then 


V/c > 1, l{k) = u{k) -L u{k — 1) -f 1. 
Then, with notation Xi = 1 — Q/N 

CTl 0-1-1 


p,(o =-E[n4‘’]=-E[n 


IW u{<Pl)+u{ai-l) + l 
•^k ‘‘'01 


k=l 

Ol-l 


k=l 


E[ n - I)]]- 


k=l 


(4.7) 


(4.8) 
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Since E[z^] = l/{2 - z), we have 




We set a{k) = 0 for fc > cti, whereas for any k < ai 


(2 


(4.9) 


SO that 


and by induction, we obtain 


;-l) _ 

2 — Xke~°-A) ’ 

(4.10) 

CTl-l 



(4.11) 

k=l 


= exp ( -^a{k)). 

(4.12) 


k>0 


Note that (4.10) reads for 1 < fc < ai 

g-a(fc) ga(fc-i) _ and + —-(4.13) 

Xk Xk ^/c+1 


and a{k — 1) > a{k) > 0 follows by induction from (4.10). Inequality (2.1) implies that 

(4.14) 




fc >0 


i>i 


i>i 


Proof of Lemma 2.2| , An explorer settling on the pile at site i, hits the i-th pile at a 
height between 1 and at. Knowing that it settles at height k, it has chance 1/Ck to settle on 
{i, k) since we are on the complete graph. We underestimate the probability of settling on 
CTj, if we only consider trajectories hitting only one of the {Ci, C 2 , • • • } before Hq. Thus, 


F(Explorer attaches pile i\ a) 


Aid 


1 

^ C/i 






(4.15) 




Fix h > 0, and write for the height occupation such that 

4‘ = a, and 4‘ = 0. 


We rewrite (4.15) in terms of the function P„ as follows 


F (Explorer attaches pile i 


1 ,P,(C'') 


^a'i’,(o 


!)■ 


(4.16) 
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For a given h < di, we now study the ratio Pg{(^)/Pg{(). As in (4.12) in t he pre vious 
paragraph, we write Pg{C^) = exp(—a(fc)), with d satisfying the relation (4.10) with 
in place of (. In other words, 


Wk > h, a{k) = a{k), and exp(—a(h — 1)) = 


2 _ e-a(/i)' 


(4.17) 


For k < h — 1, 


ga(fe-l) ^-a(k) _ _ _ ga(fc-l) ^-a{k) ^ 

Xk 


(4.18) 


We set 5k = a{k) — a{k), and from (4.18), we have that 6k > 0. In terms of 6k, (4.18) reads 
for fc < h — 1 

exp(4-i) - 1 = e-“('=-')-“('=)(exp(4) - l), (4.19) 

whereas for k = h — 1 we have 


exp{6h-i) - 1 = 


2{l-Xk) 


Xhi2 - exp(-a(h)) 


< 2. 


(4.20) 


Equality (4.20) implies (since a(.) < 2 by (4.18)) that for some constant kd 

6h-i > «d(1 - Xh) = 


(4.21) 


We deduce that 6k-i < 6k and since x {e^ — f)/x is increasing, (4.19) implies 

> exp ( — d{k — 1) — a{k)) > exp ( — a{k — 1) — a,{k)) > exp ( — 2a{k — 1)). (4.22) 


By induction on (4.22), and using (4.21), this implies that for each h < ai 

h-l 


P^iC 


h\ 


PAO 


h-l 

1 = exp (4) -1 > 4 

k=0 k=l 

h—1 h—1 

P6h-i y ^ exp ( ~ 2 

j=k 

CTl-l 




fc =0 

(h-l) 


N 


exp ( - 2 a(j)) 


3=0 


. ih-l) 

Pf^nCh ——exp 




i>i 


(4.23) 


Since Q < N, we have C/-^ < Cj: cind (4.23) implies (2.2). 


Remark 4.1 In order to obtain Corollary |2.3| consider simply the height occupation {(k+H, k 
N} and proceed along the exact same proof. This height occupation corresponds to {(dj — 
//)+, z = l,...,iV}. 
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Given A > 0, we define the early regime as the following subset of conhgurations. 


Xe{A) = {a: 


(4.24) 


Z=1 


2=1 


Dehne also k{A) to be exp(—2^4). 

Corollary 4.2 For any a G Xe{.A), and i G {1,..., 

P( Explorer attaches pile i\o') > k.{A) 

4.1 Upper Bound 


N' 


(4.25) 


Lemma 4.3 Consider diffusive deposition. Let i be a fixed site and a be a configuration 
such that ai < \/N. Then 


P( Explorer attaches to site i \ < k 


V 1)' 
N 


(4.26) 


with K = 1 + 0{N 


Remark 4.4 Comparing lower and upper bound (Lemmas 2.2 and 4.3) on attachment prob¬ 
ability we obtain a good control on this probability for conhgurations in the early regime, 
that is in d4(4) when A is small. 


Proof of Lemma f.3. If cxj = 0 then (4.26) is immediate with k = 1. If cxj > 1 the chances 


an explorer attaches to column z, in conhguration a, is smaller than if all columns distinct 
from i were set to zero. This is seen by coupling. First, for conhguration a, let a* denote 
the conhguration where we annihilate all columns distinct from i. In other words 

(<7*)^ = 0, when k ^ i, and (cr*)i = ai. 

Therefore, in cr*, the highest column is i with height ai > 1. Now, the event hit column i 
in (T® is the complement of the event hit the base hrst. Since Pg{a^) is the probability the 
explorer hits hrst the base in conhguration a®, by Lemma 2.1, we have that for some k > 0 


P(An explorer attaches to i|cr) <P(An explorer attaches to i|cr®) 


=i-nK)< 


which completes the proof. 


Lemma 4.5 Consider ballistic deposition. Consider a configuration a G then 


P( Explorer attaches to site i 



<71 + 1 

N 


(4.27) 
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Proof of Lemma 4^ To get attached to site i, the particle has to survive up to the time it 
reaches height cxj, and then at each step-down has a chance 1/A^ to fall on column i provided 
it has avoided the other columns. Thus 

F(An explorer attaches to site i\a) 


- maxcr 

n 

L h=ai + l 


^(i + 



fi V)+-') 

^ nJ\ 

N\ 

V N J 

V N J 

V N J J 




which completes the proof. 


4.2 Stochastic domination 

Both in ballistic and diffusive deposition, we have a simple upper bound on the probability 


of attaching to a given column (see Section 4.1), which depends only on its height. This, in 


turn, is used to bound the number of explorers necessary to increase the height by one unit 
in terms of a geometric random variable, for which everything can be computed explicitly. 
In other words, call ri the number of explorers needed so that column 1 reaches height 1. 
Let T 2 be the additional number of explorers needed to reach a height 2, and so on. Note 
that {ti > k} means that out of k explorers none of them has reached site 1. These times 
are used to control the height of column 1 after k explorers have been sent, 

VA;>1, VLf>l, Po{(Ti{k)> H) = P{ti + --- + th <k). 

We need to estimate the sum of the {xj} with the following general lemma. 

Lemma 4.6 Let t,T be stopping times with respeet to a filtration {4Fn\- Le-t ti := r and if 
9{n) is the time-shift by n units, define inductively 

Tn -.= TO 9{Ti H-h Tn-l). 

Let {fn, n eN} be independent random variables, which are also independent from {r„, n G 
N}. Assume that for positive integers f < f', we have 


P(r„ > e ,T > ri > P{fn > ^)P{T > O- 

Then, for any integer n and ^ > 0 

n n 

p(Y,n>i\T>i) >p(Y,P>i 


2 = 1 


2=1 


(4.28) 


(4.29) 


Similarly, if instead of (4.28) we have 

F(r„ > e ,T > e'l < P{fn > ^ > i') , 


n n 

2=1 2=1 


then 
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Proof of Lemma 4-6 We prove (4.29) by induction. The step n = 1 is obvious. Assume the 
inequality is true at step n — 1. Using that the variables are positive, 


n—1 


n—1 


T>e) = = T > e, A > e, T>^ 

K<^ i=l 


i=l 


K<t 


i=l 

n—1 


P[Tn>^-K, T>i-K 


n= 


K, T>K 


> Y,E\p{fr, > i - k)p{T >i-K 


K<i 


n—1 


""rWi=K, T>K 


n—1 


+ p(5^r,>e, T>e 
%=1 
n—1 

+ P(5^r,>e, T>i 


2 = 1 


Y,P[rn>i-K) X p[Y,n = K,T > A + p[Y,n> T>i 


K<i 


2 = 1 


2 = 1 


n—1 


P[^P + '^n> T > A . 


2=1 


(4.30) 


Now, we can exchange the role played by Ty and by ri H-in the previous argument, 

to use the induction hypothesis. Indeed, 

n—1 ^ n—1 

p( 5 ^A + f„>e, T>e) =X]p(r„ = ir) xP{Y^n>^-K, T>A + p{fn>A 

2 = 1 K=1 2 = 1 

n 

>p{J2fi>()p(T>o- 

(4.31) 


The proof of the opposite inequalities follows the same steps. 


5 Very Early Regime 


One important step in the cluster growth is to reach height log(A^). We cover this interme¬ 


diary step in the following proposition, even if it is included in Theorem 1.1 


Proposition 5.1 Consider diffusive deposition. There exist positive constants b and 7 such 
that almost surely, for N large 


(h ^ ^ 

maxailb- -— 

i<N ' logA^ 


> 7 log N. 


This proposition is concerned with what we call the very early regime, and it is based on 
comparison with urn models. 
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Two scales play an important role in this section: the time scale A^/log(iV), and the 
space scale log(A^)/loglog(iV). We therefore introduce notation 


Tn = 


N 

log{N) ’ 


and 


log(iV) 

log \og{N) ■ 


(5.1) 


The set of conhgurations a with maximal height lower than ylogiV and |cr| < /3T]\[ is called 
the very early regime and is denoted by Xvei'Jyff)- In other words, 


Xveh,/3) = w : |cr| < /3Tn, maxa^ < 7 log(iV)}, 

X 


and note that 


N 

Va e TVe(7,/5), < (maxcr,j,)|cr| </? 7 iV. 

^ x<N 

x=l 

If ta is the hitting time of set A, we show in this Section that there are constants b < jd and 
5 > 0 such that 


P{Tx„h.!if > < exp(-/V''). 


(5.2) 


Strategy of the proof. We divide time in two periods. In the hrst, of length T/v, a large 
number of columns, of order with 0 < a < 1, reach a height This is the content of 

Lemma 5.3, whose main ingredient is a coupling between diffusive deposition and random 


allocation. In the second period, we use the estimate of Corollary |2.3 to control the growth 
of these columns together with sending Poisson waves of explorers to ensure the growth of 
each column independently. 


Step 1: Reaching height Hjsf. The random allocation evolution is denoted by n h-//( n). 
Our hrst lemma deals exclusively with random allocation. 

Lemma 5.2 For a G [|, 1), and <5 < (1 — a)/2, we have almost surely, for N large enough 

\{x:p,,{Tm)>5Hm}\>N^. (5.3) 

Proof. Let X be a Poisson variable of parameter Tn/2. We have 

P{\{x : 7 ,(T^) > 6H^}\ < iV“) < P{\{x : v,{X) > SHr,}\ < N^) + P{X > T^). (5.4) 

Now, {t],,;{X),x = 1,... ,N} are independent Poisson variables of parameter l/21og(X). A 
tedious but simple computation gives 

P{vi{X) > 6 Hn) = (5.5) 


Now, 


N 

|{x : 7 x(-A) > 6Hn}\ = ^ 

X=1 


(5.6) 
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and by Bernstein’s ineqnality, for a < 1 — 26, and 


P{\{x : v,{X) > 6Hr,}\ < iV“) < exp(- ^' '' (5.7) 


Note also that from Chebychev’s exponential ineqnality 

P{X>Tn) <exp(-^^^T^ 

The statement follows. 


(5.8) 


In Lemma 3.7, we establish that diffnsive deposition, denoted 1 1 —)■ a{t) is more monopo¬ 
listic than random allocation. Thns, there is a conpling snch that with probability 1, when 
cr( 0 ) = r]{0), we have for any f > 0 


Wk<N, ai{t) > Y 


2=1 


2=1 


Assnme now that cr(TN) G and that 


L := |{a; : > 6Hn}\ > L' := \{x ; a,(T^) > 


Then, by onr conpling 


(5.9) 


> ^T/i(7Ar) > and y^n,(r/^) < -H^L. 


2=1 


2=1 


i = L' 


Then 


7 


log(iV)L' > d,(T^) > -H^L ^L'> 


2=1 


27 log log(A^) 


L. 


(5.10) 


Thus, for any a > 1 — 26, we have that V > iV". We therefore state the result as follows. 

Lemma 5.3 For a G [\, 1), and 5 < (1 — a)/2, we have almost surely, for N large enough, 
and for the diffusive deposition t h-)■ crit). 




(5.11) 


By applying Markov’s property at time Tjv 

> ('> + l)^^) < : ife(r„) > < Af") 

r 

+ sup (^Pa{'y(bT,t) € d'„(7,,3)) : a € \{x : u, > ^Hn]\ > 1V“). 

(5.12) 
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Step 2: Poisson waves. We realize diffusive deposition for times in [T/v, hT^] by a se¬ 
quence of Poisson waves, the fc-th wave made of explorers, and {X^^\ fc > 1} an i.i.d 
sequence of Poisson random variables with parameter xj^ going to inhnity with N. 

Our starting conhguration denoted satishes 


(7^°^ e AVe(7 ,/^), Aat := [x : <7^°^ > and \An\ > iV“. 

Let (7^^^ be the conhguration of diffusive deposition starting from after the k-th wave is 
sent, i.e.. 


= a 


i=i 


with |cr(^)| = |cr(°)|+ 


i=i 


Dehne now, using kd of Corollary 2.3 


.h'vE := kd exp(-7/3). 

We have for the diffusive deposition process for any A; = 1,2,..., if G /d) then 

pMt)) > Vt e W e An. (5.13) 

This immediately follows from Corollary |2.3[ 

Consider an auxiliary growth process a{t) which evolves on the sites of An U {0}, dehned 
iteratively as follows. Set a® = 0, and for i G An, set af'’ = ^Hn- Each explorer in the 
fc-th wave is attached to site i E An with probability 


pf{k) = Kne 


.(T. 




N 


(5.14) 


whereas site 0 grows by one with probability 1 — Yhi&Ki^PtiY). 
The following result is crucial. 


Lemma 5.4 There exists a coupling between t i—)■ a{t) and t i—)■ a{t) such that if < 
'Tx„e{-y,PT ^ '^dhin the k-th wave, i.e., t G [|(7 ^^“^^|,-|- |(7^^“^^|[ 

V/ G An, ai{t) > di{t). (5.15) 

Moreover, {di, i G An} are independent, where we used the shorthand notation di = 

{dY\k< 

Proof The coupling part is simple and we omit it here. We denote by {Yf, i G A^v} inde¬ 
pendent Poisson variables of parameter XNpfY)- We denote by Qi the sigma-held generated 
by XA) and by {Yf, i G An}. We now build Qk by induction. Assume that Gk-i has been 
built. Then conditioned on Qk-i, we hx the height of all sites after the k — 1-th wave. Draw 
a Poisson variable X^^^ independent of Gk-i, and denote by {Yf, i G An} the independent 
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Poisson variables of parameter pf(k)x]v which is itself Qk-i measurable. Note that de¬ 
pends only on and on the past through Y^~^, the height of site i after the 

k — 1-th wave. In other words, for any real function /j, there is a function such that 

E [/.(F‘ + --- + Yt) Ifo-i] = E [n(Yf + ... + Yt)\Y^ + ---+ if-‘] = + ... + Yp'). 

(5.16) 

Note also, that if we integrate only over X^., and for any real functions /*, for i G Ajv we 
have 

n /i(r,‘ + '" + y‘)ia*-i] = n B[/i(y‘ + --- + yflir.' + '" + Ky'] 

*SAiv isAjv /p, .. _x 

= n +■ ■ ■+ Y')- 

This means that what happens on different sites of A^r is independent. I 


Now, each Poisson wave we send has about xn explorers, and we expect to send about 
(6 — l)T;v/xAr waves. Recall that for N large, we have a.s. that G A:’ve(7,/d) and 
IAtvI > A^". Therefore, if t^ denotes the integer part of bTisf/{2exN), and for simplicity 
H=lH^ 


In tjv 

iSAjv 

k=l k=l 

In 

n (l - + '^Yt) > 'iN)). 

iGA.]^ k=l 


(5.18) 


Step 3: Dealing with one site. We show that for a function e{'y) going to 0 with 7, for 
all i G Aat 

In 

p(y^yt>l log( A^) - if) > exp (^ - e(7) log N^. (5.19) 

k=l 

We define the successive wave numbers at which the column at 1 grows. Let r be the number 
of waves needed so as to increase by at least one the height of site 1. Then, let ti = r and 
Tn = T o 9 {ti • -t- r„_i). Note that for any integer n 


F(t, > n) = F(if > = 0,..., f)'”> = 0) = exp ( - 


(5.20) 


where we used (5.14). Note that at the number of waves t = ti- 
is larger or equal than H Y k — 1. We have, using Lemma 


5.4 


'-|-Tfc_i, the configuration 


r>( \r \ ^ + k - VfxN-^ 

P[Tk > n\gk-i) < exp ( - n - — - ). 


(5.21) 


Then, we are in the setting of Lemma [T6 and have a comparison with independent geometric 
random variables {fk, k > 1} with 


E[fj,] = 


N 


K^EiH + k-iyxN 
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Then, (with the abuse of taking 7 log(A^) to be integer) 

tjv 


>7log(iV) -h r^iogiv < ^iv) > ^(^1 H-h r^iogw < ^iv) 


k=l 


7 log N 


> 


n Ln 


< 


tN 


7 log N 


k=l 
7 log N 


'jlogN 


) > (1 - exp ( - t AT 


> 


n b 


— exp ( — 


k=l 


k=l 

7&ibvE 

2e 7 ^ log(iV) 


KyE{H + kyxN 
jNlog{N) 


)) 


)) 


Now use the estimate 1 — e ^ > e “xforO<a;<a with a = 76 it'vE/ 2 e. Now, 


7 log N 


k=l 


n i-e='p( 


'ybKyE 

2 e 7 ^ log(A^) 


/ I 


( 7 log(iV ))!2 


(72 log2(Ar))7iog(iv) 

>exp(-e( 7 ) log(Ar)), 


(5.22) 


(5.23) 


where we took e( 7 ) = 7 log(ae “) — 27 . 

The proof of the Proposition 5.1 is completed if a > 6 ( 7 ). 


6 Growing Columns 


In this section we present a simple way to bound the height of the maximal pile, based on 
stochastic domination, for both ballistic and diffusive deposition. Indeed, for both models 
we construct a sequence of inter-arrival times of explorers on a given column, say column 
number 1 , stochastically dominated by independent geometric variables. 


We now state three propositions that are crucial in the proofs of Theorems 1.2 and 1.1 


The propositions bound the probabilities of building a high pile, and are proven at the end 
of this section. 


Proposition 6.1 Consider ballistic deposition. There is a constant k > 0 such that for all 
H > 2, and any site i G Gn, 


P(^ai{N) > H) < exp{—KH). 


( 6 . 1 ) 


Proposition |6.1| implies that, for a given column, N particles are not enough to reach a 
maximal height of order log{N)/K. 


Proposition 6.2 Consider diffusive deposition. For H < N^ and any site i G Gn, and 
X a positive integer 

X vr 


P7.(A') >H)< exp - -H), 


( 6 . 2 ) 


with K = 1 + 0{N ^P). 
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Proposition 6.2 implies that for a small, aN/ log(iV) particles are not enough to reach a 
maximal height of order log(A^). 

Finally, we consider diffusive deposition, with a conhguration in Xe{A), and with one 
distinguished site, say z, above height 7 log(A^). We show that as long as we do not leave 


the early regime, see equation (4.24) we have a fast growth. Let the time at which you 
exit the early regime. 

Proposition 6.3 Let XiliC he any positive eonstants with x < 1/2. Assume that there is 
a distinguished site, say i*, with ai* > 7 log(iV). Then, we have 




) < exp (- 7 (k(A)^ - |) log(]V)). (6.3) 


log(iV) 


6.1 Growing a Column in ballistic deposition 


Proof of Proposition 6.1. By lemma 4.5 


1 X k 


P{ti >k)>{l--) 


and in general 


P(Ti > k) > (l — so that E[p] > —■ 


N 


E[ti] > N. 

(6.4) 

E[n] > -. 

% 

(6.5) 


This implies that (6.1) is a large deviation event since 

H 


J2E[T,]>Nlog{H). 


i=l 


More precisely by Lemma |4.6| we have 

P[ti + ... + th < N) < P{fi + ... + < iV) 

with {fi, i = independent geometric variables of mean By the exponential 

Chebyshev’s inequality we get, for every A > 0 


H 


P{fi + ... + fH <N) < e^^Y[E[t 


—Xfi 


2 = 1 


Note that for a geometric variable X of mean 1/p 


£^[exp(—AX)] = 


p 


_ n - 


(1 -p)' 


When A is positive, exp (A) — 1 > A, and we have 


i?[exp(—AX)] < 1 


A + p‘ 
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Now, 


A 


We 


n H 

n E[exp{-Xn)] < n (1 - x + 

i=l i=l ' i=l 

choose \ = io/N so that by the asymptotic of the harmonic series 




H ^ io+H ^ 

i=l ' i=io+l 


1 ■ I 717 I 1 

1 _ . , ,lo + H 


^0 + 1 


So, we obtain 

P{fi H- \-fH<N) < exp - io[log “ l])- 

Now, if we choose io to be the integer part of aH, for some constant a, then 
(ti H- \-fH<N) < exp - aH{l - “ ^))- 


Pi 


If a is sufficiently small, say a = |, then log((l+a)/a) —1 > 0, and Lemma 6.1 


is established. 


6.2 Growing a Column in diffusive deposition 


Proof of Proposition^^ We follow the arguments of the previous proof. By using Lemma 4^ 
for any i > 1, and integer n 

P 

P{ri > n) > (1 - n—y, with k = 1 + 


By Lemma |4.6| , we have for any H, X 

P{ti + ... + th < X) < P{fi + ... Lth < X), 

with {fi, i > 1} independent geometric variables with E[fi] = N/{kP). 
Then, for every A > 0, by Chebyshev’s inequality 

^ 1 

P{n + --- + TH<X)<exp{XX-Xj2 x + myx ^ ' 


We set A = kH'^/N and we note that 

H , H 




E 


1=1 1=1 ^-W 


> H 


dx ^TT 

) 1 + 4 ’ 


( 6 . 6 ) 


(6.7) 


We conclude obtaining (6.2). 
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Proof of Proposition 6.3 Again, let {ri,r 2 ,..., } be the random number of explorers 
linked with growing a column at i* from a with an initial state with a,* = 7log(^)- By 
Corollary |4.2 we have for any integer m < X 


P„{Tk > m|ri,.. > A) < (1 - k{A) 


{'jlogN + k — 1)^ 
N 


= P{Tk > m) 


Therefore, by Lemma 4.6 


NX 


NX 


P^{J2^k>X\TN^^>X) <P{J2^k>X). 


k=l 


k=l 


By Chebyshev inequality. 




Pk 


1 , ® -^1 ^ + k-iy 

= 1 H-, with Pk := k[A) - 


e“^-(l-pfc) Pk-a 

assuming a := 1 — e~^ < pk- Thus, 

jggArj; ^ exp{a/(pfc — a)} . 

Hence, 

NX NX 

P f ^ fk> x\ < exp r - AX + ^ 


N 


k=l k=l 

.2 


Pk - a. 

We choose a = k(A)( 7 ^ log^ N)/N (and a > A/2) and X = CN/ logX, and we have 

NX 

E a , ^ 1 "' ax 

SO that 

7^ log^ N N TT 


= TTblogX, 


o 2 

P ( ^ > A') < exp ( - k( 21) + p- log n) 

k=l ^ 

7C TT 


= exp (^- 7 logX[/s:(A)^ - 


6.3 Growing a Tower in diffusive deposition 

In this section, we bound the probability of forming a high tower of explorers in diffusive 
deposition. 

Fix a region A C {!,..., X} of size L. Let H he a. fixed positive integer, and C = 
A X {0,...,X}. Given a G we define 

(T A C := {i G A : CTj < H} C Gn, and note that \a /\C\ = E (6,8) 

isA 
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Proposition 6.4 Consider diffusive deposition. For any positive H, X andf with X < N/2, 
we have that 

P(VxeA, ff,(A:)>//)<exp(L//(?“-^log(ii|±^))) (6.9) 

Moreover, for any positive a, 7 , and a positive real x satisfying 4% < 07 ^ exp(— 207 ), and the 
choice H = x^og^N), X = aT]\f and L = we have 


P{\{i : ai{aTN) > 7 logA^}| > < exp ( - xN^ ^^)- 


( 6 . 10 ) 


Proof. Lemma 4.3 immediately yields 


P( Explorer attaches to a AC \ a) < ka 


az +I 


cSctAC 


N 




Note that 

+1 < L+if 

x&af\C x&af\C 

We dehne nc{cr) = Yhx&ahC that, for N large enough we have that 


P( Explorer attaches to a A C | a) < ^[L + Hncicr)]. 


( 6 . 11 ) 


This allows us to dehne, as before, geometric random variables stochastically smaller than 
the number of explorers needed to settle one of them in C. Let ti be the number of explorers 
needed in order that one settles in C, when we start with the empty conhguration. By 
induction, when k — 1 explorers are settled in C, dehne Tk to be the number of explorers 
needed to settle the kXh explorer in C, and we do this up to time LH. Then for any 
conhguration a with nc(cr) = k — 1, for any positive integer m 


P{n > m I ct) > (1 - —f—A = P{fk > m). 


( 6 . 12 ) 


We invoke again Lemma 4^ to obtain 

P(ri H- hTHL<X) <P(ri H- f thl < X) < JJ P[exp(-Affc)] 


HL 


k=l 


HL 


< exp f AX — A y~^ 


< exp K 


-A + (/cii + P)/X 

hmx 1 


choose A = ---—f 
N 


N 




k=l 


HLi + k + Ln 


<exp(LP(^—-Oog(^^^ 7 ^ 


»)■ 


(6.13) 
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With the choice H = 'y\og{N), L = and X = aT^ we get 


P{\{i : ai{aTN) > 7 logAr}| > exp ^log( ^ ))) 

1 + e + 1/H^ 


MX 


< exp (^L log(iV/L) + LH—^ - LH^ log( ^ ^ 
<exp (^ - Llog(iV)( 7 ^ 1 og(^^) -2x-aj^^}y 


-)) 


First choose ^ = exp(— 207 ), to get 

■l + e 


(h^ log(^r^) - 2x - xh^^} > (aj^^ - 2x)} • 


Now choose 4% < ( 07 ^) exp(— 207 ), and the ineqnality (6.10) is obtained. 


(6.14) 


7 Proof of Theorems 11.11 and 11.2 


Proof of Theorem mi Proof of ( |1.9[ ): the statement follows immediately by Proposition 
with H = 31og(A^) and X < cN/ \og{N), for c small enongh. Indeed by Proposition 6.2 
have for the complementary event 


6.2 


we 


N 

P(3z : ai{X) >?,\ogN) kS" P {a i{X) > 31ogA^) 

tr (7.1) 

<N exp ^3 log N[Skc ~ ^)) • 

Hence 

P{3i : ai{X) > 3 log AT) < 

This conclndes the proof since the exponent of N is less than —1 when 9 ck < ( 37 r /4 — 2). 
(recall that k = 1 + 0 (iV“^/^)). 


Proof of (1.10). Recall that from Proposition 5.1 there is 6 > 0 (and (5.2) for the qnan- 


titative estimate), so that very likely < bT]\f, where T/v = iV/log(A^). We therefore 

condition on the evolntion np to r;t’ve( 7 ,/ 3 )‘'- 


P{ 


max a, 


faT^ + hT^) < N^) <P{rx^,{^,fsr < bTN,m&xa^{aTN + bT^) < M) 

+ -P('ra’ve(7,/3)" ^ bT^). 


(7.2) 


Now, in the hrst term nse Markov’s property at time calling for simplicity cr(T’e) : = 


P{rx,,{'y,PT < bTM,mnxax{aTN+bTN) < M) < E M 

(7.3) 
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Now, a{Xe) £ for A > 2'y[3, and there is i* G Gn such that cTi*(A’e) > 7 log(iV). 

Thus, 


max da;(aTTv) < N^) <P^(;t’e)(max(Ta;(aT 7 v) < N^\ tx,- > aTN)Pa{Xe){'rA:,‘= > o^Tn) 

+ Pa{x,){rx,<^ < ciTn, m^a^{aTN) < N^). 

(7.4) 


The first term on the right hand side of (7.4) is dealt with by Proposition 6.3 


The next lemma deals with the second term on the right hand side of ( |7.4 ). 

Lemma 7.1 Let a G Xg{A) be a configuration such that maxa^ = 7 log(iV), for some 
positive 7 . For a > 0 such that a'y < A — 1, and 4% < 07 ^ exp(— 207 ), we have 


-Po-(T'a’7 < 


maxcr. 


;(aT 7 v) < N^) < exp(— 


(7.5) 


Proof of Lemma 7.1 Note that when 07 < A — 1, 

{tx^^ < aT^, a,{aT^) < C {|{a; : a,{aT^) > 7 log(iV)}| > (7.6) 

Indeed, on the event {tx^<^ < aT^, axiaT^) < N^}, 


AN = Y1 Y1 

i i.(ri{Tx^c)<')\ogN v.ai{Tx^c)>')\ogN 

<7(logiV)aTAr + N^^\{i : ai{Tx,<^) > 7logiV}| 
^a'yN + N^^\{i : afirx.^^) > 'y\ogN}\. 


Proposition 6.4 deals with growing large towers, and by using inequality (6.10) the proof is 
complete. I 


Proof of Theorem j .4 The statement follows immediately by Proposition 6.1 with H = 
7 log(iV) with >2/k. Indeed, by Proposition 6.1, there exists a constant k such that 

P(3i : afiN) > ylogiV) < Are-^T''°g^ = 


and the proof concludes. I 

Acknowledgements. A.A. Thanks Robin Pemantle for discussions on urns. A.A. and E.S. 
thank the CIRM for a friendly atmosphere during their stay as part of a research in pairs 
program. This work has been carried out thanks to the support of A*MIDEX grant (ANR- 
ll-IDEX-0001-02) funded by the French Government ” Investissements d’Avenir” program. 
E.S. thanks Universite Paris-Est, Creteil. Finally, we thank two anonymous referees for 
their careful reviewing. 















Diffusion Limited Deposition 


34 


References 

[1] Amir, G.; One-dimensional long-range diffusion-limited aggregation III - The limit ag¬ 
gregate. Preprint (2009), arXiv:0911.0122 

[2] Amir, G.; Angel, O.; Benjamini ,L; Kozma, G.; One-dimensional long-range diffusion- 
limited aggregation I. Preprint (2009), arXiv:0910.4416 

[3] Amir, G.; Angel O.; Kozma G.; One-dimensional long-range Diffusion Limited Aggre¬ 
gation II: the transient case Preprint (2013), arXiv: 1306.4654 

[4] Asselah A., Gaudilliere A., From logarithmic to subdiffusive polynomial fluctuations for 
internal DLA and related growth models The Annals of Probability. Volume 41, Number 
3A (2013), 1115-1159. 

[5] Asselah A., Gaudilliere A., Sub-logarithmic fluctuations for internal DLA. The Annals 
of Probability. Volume 41, Number 3A (2013), 1160-1179. 

[6] Asselah, A.; Gaudilliere, A., Lower bounds on fluctuations for internal DLA. Probability 
Theory Related Fields 158 (2014), no. 1-2, 39-53. 

[7] Barlow, M.T.; Pemantle, R.; Perkins, E.A., Diffusion-limited aggregation on a tree. 
Probability Theory Related Fields 107 (1997), no. 1, 1-60. 

[8] Benjamini, I.; Duminil-Gopin, H.; Kozma, G.; Lucas, G.; The Internal Limited Aggre¬ 
gation model with random starting point. Preprint 2013. 

[9] Benjamini, L; Yadin, A., Diffusion limited aggregation on a cylinder. Gommunications 
in Mathematical Physics 279, no.1,187-223. (2008) 

[10] Diaconis, P.; Fulton, W. A growth model, a game, an algebra, Lagrange inversion, and 
characteristic classes. Rend. Sem. Mat. Univ. Politec. Torino 49 (1991), no. 1, 95-119 
(1993). 

[11] Eden,M. A two-dimensional growth process. Proc. 4th Berkeley Sympos.Math.Stat. and 
Proba, Vol IV. (1961),223-239. 

[12] Erbez-Wagner, D.; Discrete growth models. Doctorate Dissertation, University of Wash¬ 
ington, 1999. Preprint arXiv:math/9908030 

[13] T.C.Hasley, Diffusion-Limited Aggregation: A Model for Pattern Formation. Physics 
Today, 53, 11-36 (2000) 

[14] H. Kesten, How long are the arms in DLA? (1987) J. Phys. A 20, no.l, 29-33. 

[15] H. Kesten, Hitting probabilities of random walks on I/. (1987) Stochastic Processes and 
Applications 25, 165-184 

[16] H. Kesten, Upper bounds for the growth rate of DLA. (1990), Physica A 168, no.l, 
529-535. 






Diffusion Limited Deposition 


35 


[17] H. Kesten, M.V.Kozlov, F.Spitzer; A limit law for random walk in a random environ- 
mentA Compositio Mathematica, tome 30, n2 (1975), 145-168. 

[18] Jerison, D.; Levine, L.; Sheffield, S., Logarithmic fluctuations for internal DLA. 
J.Amer.Math. Soc. 25 (2012), no. 1, 271-301. 

[19] Jerison, D.; Levine, L.; Sheffield, S., Internal DLA in Higher Dimensions. 

Electron. J. Probab. 18 (2013), No. 98, 14 pp. 

[20] Jerison, D.; Levine, L.; Sheffield, S., Internal DLA and the Gaussian Free Field. 

Duke Math. J. 163 (2014), no. 2, 267-308. 

[21] Lawler, G., Intersection of Random Walks Probability and its Applications. Birkhauser 
Boston, Inc., Boston, MA, 1991. 

[22] Lawler, G.; Bramson, M.; Griffeath, D. Internal diffusion limited aggregation. Ann. 
Probab. 20 (1992), no. 4, 2117-2140. 

[23] Lawler, G; Limic, V. Random Walk: A Modern Introduction Gambridge Studies In 
Advanced Mathematics, 2010. 

[24] Levine, L.; Peres Y. Strong spherical asymptotics for rotor-router aggregation and the 
divisible sandpile. Potential Analysis 30 (2009), 1-27. 

[25] Meakin, P.; Deutch J.M. The formation of surfaces by diffusion limited annihilation 
J.Ghem.Phys. 85 (4), 1986, 2320-2325. 

[26] Norris, J.; Turner, A.; Hastings-Levitov aggregation in the small-particle limit. Gomm. 
Math. Phys. 316 (2012), no. 3, 809-841. 

[27] Pemantle, R.; A survey of random processes with reinforcement Probability Surveys, 
Vol.4 (2007), 1-79. 

[28] Richardson, D. Random growth in a tessellation. Proc. Gamb. Phil. Soc. 71 (1973), 
515-528. 

[29] Rolla, L.; Sidoravicius, V.; Absorbing-state phase transition for driven-dissipative 
stochastic dynamics on Z. Invent. Math. 188 (2012), no. 1, 127-150. 

[30] Sidoravicius, V.; Teixeira, A.; Absorbing-state transition for Stochastic Sandpiles and 
Activated Random Walks. Preprint 2014, arXiv: 1412.7098 

[31] Witten,T.A. ;Sander, L.M.; Diffusion-limited aggregation, a kinetic critical phenomenon 
(1981), Phys.Rev.Letters 47, 1400-1403. 



Diffusion Limited Deposition 


36 


A. Asselah 

Aix-Marseille Universite, Marseille, 

and LAMA, Universite Paris-Est Creteil, France. 

E.N.M. Cirillo, 

Dipartimento di Scienze di Base e Applicate per I’Ingegneria, 
Sapienza Universita di Roma, Roma, Italy. 

E. Scoppola, 

Dipartimento di Matematica e Fisica, 

Universita di Roma 3, Roma, Italy. 

B. Scoppola, 

Dipartimento di Matematica, 

Universita di Tor Vergata, Roma, Italy. 



